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I review recent applications of two-particle intensity interferometry in high energy 
physics, concentrating on relativistic heavy ion collisions. By measuring hadronic 
single-particle spectra and two-particle correlations in hadron-hadron or heavy-ion 
collisions, the size and dynamical state of the collision fireball at freeze-out can 
be reconstructed. I discuss the relevant theoretical methods and their limitations. 
By applying the formalism to recent pion correlation data from Pb+Pb collisions 
at CERN we demonstrate that the collision zone has undergone strong transverse 
growth before freeze-out (by a factor 2 in each direction), and that it expands 
both longitudinally and transversally. From the thermal and flow energy density at 
freeze-out the energy density at the onset of transverse expansion can be estimated 
from conservation laws. It comfortably exceeds the critical value for the transition 
to color deconfined matter. 



1 Introduction 

The principle of two-particle intensity interferometry, developed by R. Han- 
bury Brown and R.Q. Twiss in the mid 1950's for radio and optical astronomy i 
was independently rediscovered by particle physicists just a few years latere 
It has since become known under the names "HBT interferometry" , "GGLP 
effect" , and "Bose-Einstein correlations" . The method exploits the effects on 
the phase-space density of Bose-Einstein symmetrization (or Pauli antisym- 
metrization) of multiparticle states of identical particles. In astronomy one 
studies photon pairs, in nuclear and particle physics one typically uses pion, 
kaon or nucleon pairs (recently also pion tripletsHu) Strictly speaking, the 
method is used oppositely in astronomy and high energy physics: in astronomy 
one measures the two-photon correlation function as a function of the space- 
time distance of the two detected photons and extracts from it information 
about the size of the emitter in momentum space (specifically in the opening 
angle between the two photon momentum vectors). By knowing the distance 
of the emitter (star, radio source) this "angular size" in momentum space can 
then be translated into a spatial radius of the source. In particle physics one 
measures the correlation function as a function of the momentum difference 
between the two pions and extracts from it information about the space-time 
extension of the emitting source. These two opposite ways of looking at inten- 
sity interferometry emphasize its generic nature as a phase-space effect which 
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can only be understood if both the space-time and momentum-space structure 
of the emitter are taken into account simultaneously. 

Another big difference between HBT interferometry in astronomy and nu- 
clear physics is that the sources studied in astronomy are static on the timescale 
of the observation, while the sources created in high energy collisions are highly 
dynamical and very shortlived. As I will show this has dramatic implications 
for the HBT formalism to be used for the analysis of measured correlation 
functions which were only gradually realized during the last 10 years and for 
which a common understanding was developed only quite recently. 

Relativistic heavy ion collisions are performed in order to create color de- 
confined strongly interacting matter, the "quark-gluon plasma" . This plasma 
is hot and generates a huge pressure which drives a strong collective expan- 
sion of the reaction zone into the surrounding vacuum. As a consequence of 
expansion, the quark-gluon matter cools and dilutes until it can no longer 
remain in the color deconfined state: it hadronizes. If the reaction zone is suf- 
ficiently large and locally equilibrated, this process manifests itself as a bulk 
phase transition, very similar to the quark-hadron phase transition in the very 
early universe (about 20-40 fia after the Big Bang). Two crucial questions are 
therefore (a) whether the initial energy density in the reaction region is large 
enough to make color deconfined quark-gluon matter, and (b) whether this 
partonic system equilibrates sufficiently to undergo a thermodynamic phase 
transition while expanding. While the longitudinal expansion of the reaction 
zone may be simply due to incomplete stopping of the two colliding nuclei, 
transverse collective expansion flow can only be generated by the buildup of a 
locally isotropic pressure component. Since this requires a certain degree of lo- 
cal equilibration of the momentum distributions, the identification of collective 
transverse flow, superimposed by random thermal motion, plays an essential 
role in any attempt to answer these two questions. As we will see, two-particle 
interferometry features prominently in such an endeavour. 

Over the last years ample evidence was accumulated that the hot and 
dense collision region in relativistic heavy ion collisions indeed thermalizes and 
shows collective dynamical behaviour. Most of this evidence is based on a 
comprehensive analysis of the hadronic single particle spectra. It was shown 
that all available data on hadron production in heavy ion collisions at the 
AGS and the SPS can be understood within a simple model which assumes 
locally thermalized momentum distributions at freeze-out, superimposed by 
collective Jxyrlrodynamical expansion in both the longitudinal and transverse 
directionsHtj The collective dynamical behaviour in the transverse direction is 
reflected by a characteristic dependence of the inverse slope parameters of the 
m^-spectra ("effective temperatures") at small mj_ on the hadron massesu 
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New data from the Au+Au and Pb+Pb systems support this picture and 
show that the transverse collective dynamics is much more strongly exhibited 
in larger collision systems than in the smaller ones from the first rounds of 
experiments. The amount of transverse flow also appears to increase mono- 
tonically with collision energy from GSI/SIS to AGS energies, but may show 
signs of saturation at the even higher SPS energyB 

The extraction of flow velocities and thermal freeze-out temperatures from 
the measured single particle spectra relies heavily on model assumptions!! The 
single-particle spectra are ambiguous because they contain no direct informa- 
tion on the space-time structure and the space-momentum correlations induced 
by collective flow. In terms of the phase-space density at freeze-out ("emis- 
sion function") S(x,p) the single-particle spectrum is given by EdN/d 3 p = 
J d 4 x S(x,p); the space-time information in S is completely washed out by 
integration. Thus, on the single-particle level, comprehensive model studies 
are required to show that a simple hydrodynamical model with only a few 
thermodynamic and collective parameters can fit all the data, and additional 
consistency checks are needed to show that the extracted fit parameter values 
lead to an internally consistent theoretical picture. The published literature 
abounds with examples demonstrating that without such consistency checks 
the theoretical ambiguity of the single particle spectra is nearly infinite. 

At this point Bose- Einstein correlations between the momenta of identical 
particle pairs provide crucial new input. They give direct access to the space- 
time structure of the source and its collective dynamics. In spite of some 
remaining model dependence, the set of possible model sources can thus be 
reduced dramatically. The two-particle correlation function C(q, K) is usually 
well approximated by a Gaussian in the relative momentum q whose width 
parameters are called "HBT (Hanbury Brown-Twiss) radii". It was recently 
shown that these radius parameters measure certain combinations of 

the second central space-time moments of the source. In general Jjhey mix 
the spatial and temporal structure of the source in a nontrivial way£3 and the 
remaining model dependence enters when trying to unfold these aspects. 

Collective dynamics of the source leads to a dependence of the HBT radii 
on the pair momentum ^P F^ Th is feature was recently quantitatively reana- 
lyzed, both analyticallyB0llM3 and numericallyElllZf The velocity gradients 
associated with collective expansion lead to a dynamical decoupling of different 
source regions, and the HBT radii measure the size of the resulting "space-time 
regions of homogeneity" of the sourcelH'E[H around the point of maximum emis- 
sivity for particles with the measured momentum K. The velocity gradients 
are smeared out by random thermal motion of the emitters around the fluid 
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velocity^ Due to the exponential decrease of the Maxwell distribution, this 
thermal smearing factor shrinks with increasing transverse pair momentum 
K±; this is the basic reason for the if^-dependence of the HBT radii. 

Unfortunately, other gradients in the source (for example spatial and tem- 
poral temperature gradients) can also generate a if-dependence of the HBT 
radii! If ^| Furthermore, the pion spectra in particular are affected by res- 



onance decay contributions, but on ly at s mall K±. This may also affect the 
HBT radii in a Xj^-dependent wayllMJJa The isolation of collective flow, in 
particular transverse flow, from the i^-dependence of the HBT radii thus 
requires a careful investigation of these different effects. 

Our group studied this if-dependence of the HBT radii within a simple 
analytical model for a finite thermalized source which expands both longitudi- 
nally and transversally. For presentation I use the Yano-Koonin-Podgoretskii 
(YKP) parametrization of the correlator which, for sources with dominant 
longitudinal expansion, provides an optimal separation of the spatial and tem- 
poral aspects of the sourcet3cZI The YKP radius parameters are independent 
of the longitudinal velocity of the observer frame. Furthermore, in all thermal 
models without transverse collective flow, they show perfect Mj_-scaling (in 
the absence of resonance decay contributions) . Only the transverse gradients 
induced by a non-zero transverse flow can break this Mj_-scaling, causing an 
explicit dependence on the particle rest mass. This allows for a rather model- 
independent identification of transverse flow from accurate measurements of 
the YKP correlation radii for pions and kaons. High-quality data should also 
allow to control the effects from resonance decays. 

A comprehensive and didactical discussion of the formalism and a more 
extensive selection of numerical examples can be found in the lecture notes EH 
to which I refer the interested reader for more details. 

2 Spectra and emission function 

2.1 Single-particle spectra and two-particle correlations 

The covariant 1- and 2-particle momentum spectra arc defined by 



dN 
d 3 p 

dN 



= E ~ = E(a+a p ), (1) 



Pi (Pa ,Pb) = E a E b ^ = E a E b (d+ a o+ d Pb a Pa ) , (2) 

where a+ (a p ) creates (destroys) a particle with momentum p. They are 
normalized to (N) and (N(N — 1)) (i.e. the average number of pions or of pions 
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in pairs per event), respectively. The angular brackets denote an ensemble 
average (O) = tr (pO) where p is the density operator associated with the 
ensemble. The two-particle correlation function is defined as 

If the two particles arc emitted independently and final state interactions are 



on 

neglected one can provelfHI a generalized Wick theorem 



C{Pa,Pb) = 1± + — + ; ■ (4) 

Note that the second term is positive definite, i.e. the correlation function 
cannot, for example, oscillate around unity. This is no longer true if final 
state interactions are included (see below) . I will here assume that the emitted 
particles are bosons which I will call pions. 



2.2 Source Wigner function and spectra 

In the language of the covariant current formalismEU the source of the emitted 
pions can be described in terms of classical currents J(x) which act as classical 
sources of freely propagating pions. They parametrize the last collision from 
which the free outgoing pion emerges. Very helpful for the following will be 
the so-called "emission function" S(x. K): t-jlc3 

S(x,K) = J JLV e - iK -y (r{x+±y)J{x-\y)) . (5) 

It is the Wigner transform of the density matrix associated with the classical 
source amplitudes J(x). This Wigner density is a quantum mechanical object 
defined in phase-space (x, K); it is real but not always positive definite. Text- 
books on Wigner functions show that their non-positivity is a genuine quan- 
tum effect resulting from the uncertainty relation and is concentrated at short 
phase-space distances; when the Wigner function is averaged over phase-space 
volumes large compared to the volume (27r?i) 3 of an elementary phase-space 
cell, the result is real and positive definite and behaves exactly like a classical 
phase-space density. 

The emission function S(x, K) is thus the quantum mechanical analogue 
of the classical phase-space distribution which gives the probability of finding 
at point x a source which emits free pions with momentum K. It allows to 
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express the single-particle spectra a pd tw o-particle correlation function via the 
following fundamental relations:!! 

E p^T = fd 4 xS(x,p), (6) 



d 3 p 

\J d 4 xS(x,K)e l ^ x \ 2 
C(q,K) - 1 + j — — R+ i q) j dix ^ K-\q)' (?) 

For the single-particle spectrum (|^), the Wigner function S(x,p) on the r.h.s. 
must be evaluated on-shell, i.e. at p° = E p = y/ m 2 + p 2 . For the correlator (0) 
we have defined the relative momentum q — p a —pbi Q = E a — Eb between the 
two particles in the pair, and the total momentum of the pair K = (p a +pb)/2, 
K° = (E a + E},)/2. Of course, since the 4-momenta p a _b of the two measured 
particles are on-shell, = Ei = y m 2 + p 2 , the 4-momenta q and K are in 
general off-shell. They satisfy the orthogonality relation 

q-K = m 2 a -m 2 b =0. (8) 

Thus, the Wigner function on the r.h.s. of Eq. (^) is not evaluated at the 
on-shell point K° = Ejc- This implies that for the correlator, in principle, 
we need to know the off-shell behaviour of the emission function, i.e. the 
quantum mechanical structure of the source! Fortunately, nature is nice to 
us: the interesting behaviour of the correlator (its deviation from unity) is 
concentrated at small values of |q|. Expanding K° = (E a + Eb)/2 for small q 
one finds 

Since the relevant range of q is given by the inverse "size" of the source (to be 
defined more exactly below) , this approximation is valid as long as the Comp- 
ton wavelength of the particles is small compared to this "source size". For 
pion, kaon, or proton interferometry in heavy-ion collisions this is the case, due 
to the particle rest masses. This is of enormous practical importance: it allows 
one to replace the source Wigner density by a classical phase-space distribution 
function for on-shell particles. This provides a necessary theoretical foundation 
for the calculation of HBT correlations from classical hydrodynamic or kinetic 
(e.g. cascade or molecular dynamics) simulations of the collision. 

If one approximates the product of single-particle distributions in the de- 
nominator of (0) by the square of the spectrum at the average momentum K, 
one obtains the simpler result 



J d^xe 1 ^ S(x,K) 



J d 4 x S(x, K) 



= l + \(e lq - x )\ . (10) 
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The deviations from this approximation are proportional to_the curvature of 
the single-particle distribution in logarithmic representation^ They are small 
in practice because the measured single-particle spectra are usually more or 
less exponential. In the second equality of (|l^) we defined (. . .) as the average 
taken with the emission function; due to the independence of S(x, K) this 
average is a function of K. This notation will be used extensively below. 

The fundamental relations (|^J^,|l^) show that both the single-particle spec- 
trum and the two-particle correlation function can be expressed as simple in- 
tegrals over the emission function. The emission function thus is the basic 
ingredient in the theory of HBT interferometry: once it is known, the calcula- 
tion of one- and two-particle spectra is straightforward (even if the evaluation 
of the integrals may in some cases be technically involved). More interest- 
ingly, measurements of the one- and two-particle spectra provide access to the 
emission function and thus to the space-time structure of the source. This 
latter aspect is, of course, the motivation for exploiting HBT in practice. In 
my talk I will concentrate on the question to what extent this access to the 
space-time structure from only momentum-space data really works, whether it 
is complete, and (since we will find it is not and HBT analyses will thus neces- 
sarily be model-dependent) what can be reliably said about the extension and 
dynamical space-time structure of the source anyhow, based on a minimal set 
of intuitive and highly suggestive model assumptions. 

2.3 Final state interactions (FSI) 

Equation (]?]), via the plane wave factor e lq ' x under the integral of the exchange 
term, reflects the absence of final state interactions: 

P2(Pa,Pb)= f [5(xi,iif+§) S(x 2 , K-%)±e iq < Xl - x ^S(x 1 , K)S(x 2 ,K)}. 

J xi ,x 2 

In practice particle interferometry is done with charged particle pairs which 
suffer long-range Coulomb final state repulsion on their way out to the detector. 
In addition, there may be strong final state interactions, e.g. in proton-proton 
interferometry where there is a strong s-wave resonance jusL_above the two- 
particle threshold. In this case Eq. (ftW must be replaced byt 
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P2(p a ,Pb)= S(x+ S(x- %,p b ) 

J x,y 

x %°) K/2(y-^ )| 2 + 0{-y°) \4> g/2 (y-v a y°)\ 2 
± I S(x + *,K) S(x-*K) <f>*_ q/2 (y-vy°) <f> q/2 (y-vy°). (11) 
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Here v=K/Ek, v a =p a /Ex, Vb=Pb/Ex are (to quadratic accuracy in q) the 
velocities of the particles with momentum K, p a , pt, respectively, and 4> q /2( r ) 
is an FSI distorted wave with asymptotic relative momentum q/2, evaluated 
at the two-particle relative distance r. Upon replacing the latter by plane 
waves ^TTJ) turns into (0). The FSI distorted waves can be calculated by 
solving a non-relativistic Schrodinger equation for the relative motion which 
includes the FSI potential in the rest system of the pair (where K = v = 0). 
Eq. ( |TT| ) represents a non-relativistic Galilei-transformation of the result from 
the pair rest frame to the frame in which p a and p\> are measured; therefore it 
is only valid in observer frames in which the pair moves non-relativistically. In 
order to evaluate Eq. ( |ll| ) one must therefore first transform the 4-momenta 
Pa.b to such a frame (best directly into the pair rest frame). The momentum 
argument q of the FSI distorted waves </> is then the difference between the two 
spatial momenta in that frame, and their space-time argument y — Viy° is the 
relative distance of t he two particles in that frame at the time when the second 
particle is emitted^S Since the latter depends not only on the time difference y° 
between emission points, but also on the velocity of the first emitted particle, 
these arguments depend on the momentum argument of the emission function 
associated with the first emitted particle. The two terms ~ 0(±y°) in the direct 
term reflect the two possible time orderings between the emission points. 

Equations (O) and ([ll]) can h e imp lemented into event generators, follow- 
ing the procedure given elsewherelfj 



2.4 The redundance of wavepackets 

In the last two years it has been repeatedly suggested that, due to the smallness 
of the sources, the theory of particle interferometry in high energy physics 
should be based on a finite-size wave-packet descripti on rather th an on plane 
wave propagation from the source to the detectorr^r^lr^r^lr^l I will argue 
that, if done correctly, it doesn't matter whether you start from wave packets 
or not, and that the free parameter a (the wave packet width) occurring in 
these approaches is essentially unmeasurable (except, perhaps, in elementary 
e + e~ or pp collisions). 

To the extent that the detector really measures the momenta of the par- 
ticles (which it is supposed to do with highest possible accuracy), the mea- 
surement process can be described as a projection at t = 00 of the emitted 
pion state on a plane wave momentum eigenstate, irrespective of the actual 
localization properties of the emitted states; the latter will be reflected in the 
momentum spectrum resulting from this projection. The usefulness of HBT in- 
terferometry rests exactly on the fact that free propagation after emission from 
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the source does not change this overlap integral, i.e. that the 1- and 2-particle 
momentum spectra remain unchanged on their way from source to detector. 
Otherwise data measured meters away from the collision could not be used 
to extract information about the reaction zone. Final state interactions spoil 
this feature; therefore they must be accounted for analytically or numerically 
via Eq. (|ll]) before source information can be extracted from the measured 
correlator. For free propagation, however, the spectrum calculated from the 
overlap between the emitted wavefunction and the momentum eigenstates is 
identical if calculated at the time of particle freeze-out or at the detector time 
several nanoseconds laterlH In both cases, the relation between the spectra 
and the source distribution S(x, K) is given by Eqs. (||) and (0). 

While wave packets are a physically intuitive concept, one must take care 
not to exaggerate their role in particle interferometry. I will show that, to first 
order, their intrinsic structure can be completely absorbed into the emission 
function, le aadgg n o measurable trace of the wave packet width. Statements to 



the contrarylfUE^I are at best misleading. To prove my point, let me anticipate 
from Sec. 3.1 that the primary feature of the two-particle correlator, its Gaus- 
sian width in q, determines only the r.m.s. width of the source in space-time; 
the extraction of finer structures of the source requires a quantitative study 
of the deviations from Gaussian behaviour of the correlator and therefore at 
least an order of magnitude more accurate data. The dominant features of 
the correlator can thus be reproduced by replacing the true source function 
S(x, K) by a Gaussian with the same r.m.s. width in space-time. 

Let us now follow custom HJU'lZHDIl and assume that the source emits 
at times Tj Gaussian wave packets of spatial width a, centered at points £i in 
coordinate space and 7Tj in momentum space: 

= (7rfj2)3/4 exp ^ + 4 7T 4 • xj . (12) 

The Wigner density corresponding to such an individual wave packet is given 
by 

S w . p .(x,p) = ^6( X °-T i ) exp (- (X ~J l)2 -<J 2 (P-^ ; (13) 

it is normalized to j(d 3 p/E p ) J d x S w , p .(x,p) = 1. Its r.m.s. width parame- 
ters Ax = <r/\/2 and Ap — l/(^/2a in the three Cartesian directions saturate 
the uncertainty relation AxAp > h/2. 

Let us distribute the wave packet emission points (Ti,£i,iVi) according to 
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the classical Gaussian phase-space distribution 



P H,)=(iv)^M^4^ (14) 

normalized to the average pion multiplicity (N) per event. The r.m.s. widths 
At, Ro, Aq of this classical probability density are not constrained by quan- 
tum mechanics. It can be shown 0^301^1 that the spectra and correlation 
functions are then given by Eqs. (0) and (0) with an emission function S(x, K) 
which is obtained by folding the classical phase-space distribution (|l4| ) with 
the intrinsic Wigner density of the wave packets: 

S(x, K) — J dr d 3 C d 3 n P (t, £, tt) 5 w . p . (x° — t,x - K - tv) 

exp I 2 (At) 2 J ex P 1 2JF 2A^ J . s 

- w OT^ U ■ (15) 



R 2 = Ris + y, A a = A2 + ^. (16) 

Only the combinations R, A from ( p^ ) enter the 1- and 2-particle spectra; at 
the Gaussian level there exists no measurement which allows to disentangle 
Ro, A and a. As long as RA 3> h/2, there are infinitely many combinations 
of Ro, Ao,er which describe the same source. Only if R and A saturate the 
uncertainty limit, RA — Ti/2, one may conclude Ro = A = and a = \/2R. 

The wave packet size a is thus generally not measurable; as long as the 
effective emission function S(x, K) (i.e. R, A) is not changed, variations of a 
have no influence on the momentum distributions. Final state interactions do 
not affect this reasoning either - instead of (0) one must then evaluate ( [ll|) with 
the same emission function (|l5|). The effect of 2-particle final state interactions 
on the time evolution of the wave packets c3 is fully taken into account by 
weighting the effective emission function S(x, K) with the distorted waves 4> q /2 
in (|TT|) ; a generalization which includes also the FSI with the electric charge 
of the central fireball is known 0111. There is no need to describe the time 
evolution of the wave packets explicitly as done by Merlitz and PelteEfl using 
a sophisticated numerical algorithm. Finally, the above statement remains 



1 



true 



32 



are included. 



even if multi-boson symmetrization effects t 
I have presented the argument using Gaussian wave packets and a Gaussian 
parametrization for the distribution of their centers. As explained above this 
is sufficient since, in leading order, HBT interferometry gives only access to the 
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r.m.s. widths of the source such that one cannot distinguish Gaussians from 
other source shapes. At the next level of accuracy these statements are surely 
subject to (small) corrections. So far, however, I cannot see how in heavy-ion 
collisions the wave packet size might be determined experimentally. For this 
reason I prefer avoiding the concept of wave packets, thereby eliminating the 
poorly controlled free parameter <j from the theory of HBT interferometry in 
nuclear physics. 

From here on I will neglect final state interactions, assuming that the 
measured correlators have been or can be corrected for them. 



2.5 The mass-shell constraint 



Expressions ([7|jlC|) show that the correlation function is related to the emission 
function by a Fourier transformation. At first sight this might suggest that 
one should easily be able to reconstruct the emission function from the mea- 
sured correlation function by inverse Fourier transformation, the single particle 
spectrum (^|) providing the normalization. This is, however, not correct. The 
reason is that, since the correlation function is constructed from the on-shell 
momenta of the measured particle pairs, not all four components of the relative 
momentum q occurring on the r.h.s. of ( |Io|) are independent. They are related 
by the "mass-shell constraint" (|J) which can, for instance, be solved for q°: 

q °=f3.q with 0= * (17) 

j3 is (approximately) the velocity of the cm. of the particle pair. The Fourier 
transform in ( |lC)| ) is therefore not invertible, and the reconstruction of the 
space-time structure of the source from HBT measurements will thus always 
require additional model assumptions. 

It is instructive to insert ( |l7|) into (|l0|): 



J d 4 x exp(iq-(x- (3t)) S{x,K) 



Jd A xS(x,K) 



(18) 



This shows that the correlator C(q, K) actually mixes the spatial and temporal 
information on the source in a non-trivial way which depends on the pair 
velocity (3. Only for a pulsed source ~ 6(t—to) things are simple: the correlator 
then just measures the Fourier transform of the spatial source distribution. 

It is instructive to look at the problem also in the following way: If one 
rewrites Eq. (n8h in the pair rest frame where K — and hence q° — 0, one 
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obtains 



where 



with 



C(q,K)-l = J d\ cos{q-r)S Ie i{r;K) (19) 



S iel (r;K) = J d 3 Rs K (R+±r) s K (R~ \r) , 



(20) 



is the time-integrated normalized relative distance distribution in the source. 
The la tter can, in principle, be uniquely reconstructed from the measured cor- 



relator 153 by inverting the cosine- Fourier transform ([191). But since it gives 



only the time integral of the relative distance distribution for fixed pair mo- 
mentum K in the pair rest frame, no direct information on the time structure 
of the source is obtainable! Only by looking at the result as a function of 
K, which, as I will show, brings out the collective dynamical features of the 
source, can one hope to unfold the time-dependence of the emission function. 
It is clear that this will be only possible within the context of specific source 
parametrizations . 



2.6 K-dependence of the correlator 

We have seen that in general the correlator is a function of both q and K. 
Only if the emission function factorizes in x and K, S(x,K) = F{x) G(K) 
(no "x-if -correlations" , i.e. every point x in the source emits particles with 
the same momentum spectrum G(K)), the if-dependence in G(K) cancels 
between numerator and denominator of ([io|), and the correlator seems to be K- 
independent. However, not even this is really true: even after the cancellation 
of the explicit ^-dependence G{K), there remains an implicit if-dependence 
via the pair velocity (3 « K/E^ in the exponent on the r.h.s. of Eq. (plf)! 
Only if both conditions, factorization of the emission function in x and K 
and instantaneous emission ~ 5(t — t a ), apply simultaneously, the correlation 
function is truely if -independent. Even for stars which are time- independent 
sources there remains a if-dependence: the correlator is thenl§3 ~ S((3 ■ q), i.e. 
there are only transverse (_L (3), but no longitudinal (|| (3) correlations. 

It is hard to believe that this complication in the application of the original 
HBT idea to high-energy collisions went nearly unnoticepLfor more than 20 
years. It was only brought to light in 1984 by Scott Pratt 1111 in his pioneering 
work on HBT interferometry for heavy- ion collisions. 
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If one parametrises the correlator by a Gaussian in q (see below) this means 
that in general the parameters ("HBT radii") depend on K. Typical sources 
of x-K correlations in the emission function are a collective expansion of the 
emitter and/or temperature gradients in the particle source: in both cases 
the momentum spectrum ~ exp[— p-u(x)/T(x)] of the emitted particles (where 
u^{x) is the 4- velocity of the expansion flow) depends on the emission point. In 
the case of collective expansion, the spectra from different emission points are 
Doppler shifted relative to each other. If there are temperature gradients, e.g. 
a high temperature in the center and cooler matter at the edges, the source 
will look smaller for high-momentum particles (which come mostly from the 
hot center) than for low-momentum ones (which receive larger contributions 
also from the cooler outward regions). 

We thus see that collective expansion of the source induces a if-dependence 
of the correlation function. But so do temperature gradients. The crucial ques- 
tion is: does a careful measurement of the correlation function, in particular of 
its .ftf-dependence, permit a separation of such effects, i.e. can the collective dy- 
namics of the source be quantitatively determined through HBT experiments? 
We will see that this is not an easy task; however, with sufficiently good data, 
it should be possible. In any case, the -fT-dependence of the correlator is a 
decisive feature which puts the HBT game into a completely new ball park. 
Two-particle correlation measurements which are not able to resolve the in- 
dependence of the HBT parameters arc, in high energy nuclear and particle 
physics, only of very limited use. 

3 Model-independent discussion of HBT correlation functions 

3.1 The Gaussian approximation - HBT radii as homogeneity lengths 

The most interesting feature of the two-particle correlation function is its half- 
width. Actually, since the relative momentum q = p\ — p2 has three Cartesian 
components, the fall-off of the correlator for increasing q is_not described by a 
single half- width, but rather by a (symmetric) 3x3 tensor 1121 which describes 
the curvature of the correlation function near q = 0. We will see that in 
fact nearly all relevant information that can be extracted from the correlation 
function resides in the 6 independent components of this tensor. This in turn 
implies that in order to compute the correlation function C it is sufficient to 
approximate the source function S by a Gaussian in x which contains only 
information on its space-time moments up to second order. 

Let us write the arbitrary emission function S{x, K) in the following form: 

S(x, K) = N(K) S(x(K),K) e~^ {K) B ^ {K) *" (K) + SS(x, K) , (22) 
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where we adjust the parameters N(K), x^(K) 1 and B^ V (K) of the Gaussian 
first term in such a way that the correction term 6S has vanishing zeroth, first 
and second order space-time moments: 



fx 8S(x, K) = J (fx x" 5S{x, K) = J (fx x"/ 8S(x, K) = . (23) 
This is achieved by setting 

dN del 1 ' 2 B, V (K) , i 

N{K) - EK ^K S(x(K),K) ' (24) 

x»(K) = (x"), (25) 

{B-^^JK) = (x^^dx-xUx-x),). (26) 

The (if-dependent) average over the source function (. . .) has been defined 
in Eq. (|l0|). The normalization factor (Q) ensures that the Gaussian term in 
(^2|) gives the correct single-particle spectrum (^) ; it fixes the normalization on- 
shell, i.e. for K° = Ek, but as we discussed this is where we need the emission 
function also for the computation of the correlator. x(K) in ( pBj ) is the centre 
of the emission function S(x, K) and approximately equal to its "saddle point" , 
i.e. the point of highest emissivity for particles with momentum K. The second 
equality in ( p6|) defines x as the space-time coordinate relative to the centre of 
the emission function; only this quantity enters the further discussion, since, 
due to the invariance of the momentum spectra under arbitrary translations of 
the source in coordinate space, the absolute position of the emission point is 
not measurable in experiments which determine only partickL.momenta. Since 
x(K) is not measurable, neither is the normalization A^(i^)H3 as its definition 
( pi| ) involves the emission function at x(K). Finally, Eq. ( ptf ) ensures that 
the Gaussian first term in (|2^) correctly reproduces the second central space- 
time moments {x^x^) of the original emission function, in particular its r.m.s. 
widths in the various space-time directions. 

Inserting the decomposition (^2|) into Eq. ( ]To| ) we obtain 

C(q, K) = 1 + exp [-ftf (x»x v ){K)] + SC(q, K) . (27) 

The Gaussian in q results from the Fourier transform of the Gaussian contri- 
bution in (^2~|); the last term SC receives contributions from the second term 
SS in (|22| ) which contains information on the third and higher order space- 
time moments of the emission function, like sharp edges, wiggles, secondary 
peaks, or non-Gaussian tails in the source. It is at least of order q 4 ; the second 
derivative of the full correlator at q = is given exactly by the Gaussian in 

(§)■ 
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In the past it has repeatedly been observed that the correlation data ap- 
pear to be better fit by exponentials than by Gaussians. As far as I know, 
however, this happened usually for 1-dimensional fits as a function of the sin- 
gle Lorentz invariant variable Qf nv = (q ) 2 — q 2 while, at least for heavy ion 
collisions, the 3-dimensional correlators look much more Gaussian. (Correla- 
tors from elementary collisions are more strongly affected by resonance decay 
contributions (Sec.^), givin g r ise to severe deviations from Gaussian behavior 
even in three g-dimensionsE^I) Contemplating the structure of Eq. ( |27| ) one 
realizes that a fit in Qf nv does not make sense: the generic structure of the 
exponent, —q^q u {x^x v ), tells us that the term (g ) 2 should come with the time 
variance of the source while the spatial components (q 1 ) 2 should come with the 
spatial variances of the source. Since all variances are positive semidefinite by 
definition, it does not make sense to parametrize the correlation function by a 
variable in which (q ) 2 and q 2 appear with the opposite sign! Such a fit might 
work if the time variance and all mixed variances vanished identically and all 
three spatial variances were equal, but this is certainly not generic and also not 
frame- independent. The variable Qi nv should therefore not be used for fitting 
HBT data. 

Note that Eq. ( |27| ) has no factor \ in the exponent. If the measured 
correlator is fitted by a Gaussian as defined in (|27]) , its q- width can be directly 
interpreted in terms of the r.m.s. widths of the source in coordinate space. 
Model comparisons are thus most easy if the latter are directly parametrized 
in terms of r.m.s. widths. 

Eqs. and ( ^7|) would not be useful if the contributions from 5S and 
5C were not somehow small enough to be neglected. This requires a numer- 
ical investigation. It was shown EL3 that in typical (and even in some not so 
typical) situations SS has a negligible influence on the half width of the corre- 
lation function. It contributes only weak, essentially immeasurable structures 
in C(q, K) at large values of q. The reader can easily verify this analytically 
for an emission function with a sharp box profile; the results for the exact 
correlator and the one resulting from the Gaussian approximation (22) differ 
by less than 5% in the half widths^ the exact correlator has, as a function of 
q, secondary maxima with an amplitude below 5% of the value of the correla- 
tor at q = 0. We have checked that similar statements remain even true for a 
source with a doughnut structure, i.e. with a hole in the middle, which was ob- 
tained by rotating the superposition of two 1-dimensional Gaussians separated 
by twice their r.m.s. widths around their center. The only situation where 
these statements require qualification is if the correlator receives contributions 
from the decay of long-lived resonances; unfortunately, this is of relevance for 
pion interferometry as will be discussed in Sec. a 
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Eq. ( p7| ) implies that the two-particle correlator measures the second cen- 
tral space-time moments of the emission function. That's it - finer features 
of its space-time structure (edges, wiggles, holes) cannot be measured with 
two-particle correlations, but require the analysis of three-, four-, . . . , many- 
particle correlations^ It follows that possible rapid quantum oscillations of the 
source Wigner density are also essentially unmeasurable by 2-particle inter- 
ferometry. The variances (x^x^) are in general not identical with our naive 
intuitive notion of the "source radius" : unless the source is stationary and has 
no x-i"T-correlations at all, the variances depend on the pair momentum K and 
cannot be interpreted in terms of simple overall source geometry. Their correct 
interpretation BEHO is in terms of "lengths of homogeneity" which give, for 
each pair momentum K, the size of the region around the point of maximal 
emissivity x(K) over which the emission function is sufficiently homogeneous 
to contribute to the correlator. Thus HBT measures "regions of homogeneity" 
in the source and their variation with the momentum of the particle pairs. As 
we will see, the latter is the key to their physical interpretation. 



3.2 YKP parametrization for the correlator and HBT radius parameters 

A full characterization of the source in terms of its second order space-time 
variances requires knowledge of the 10 functions {x^x v ){K). However, due to 
the mass-shell constraint (ttjh which leaves only three independent components 



of q, only 6 linear combinations of theses functions are actually measurable!!^ 
For azimuthally symmetric sources three of these 10 functions vanish by sym- 
metryll^l, but again the mass-shell constraint permits to measure only 4 linear 
combinations of the remaining 7 functions of K . 

Before the correlator (^7|) can thus be fit to data, the redundant compo- 
nents of q must first be eliminated via (|l7j). We use a Cartesian coordinate 
system with the z-axis along the beam direction and the x-axis along K±. 
Then (3 = (/3i,0, /?;). We assume an azimuthally symmetric source (impact 
parameter w 0) and eliminate from ( |27j ) q x and g v in terms of q\ = q\ +q%, qi 
and q°. This yields the YKP parametrizatior@li3 



C(q, K) = 1 + exp 



-R\ql-R\ (qf-(q f)-(Rl+Rf)(q-U) 



(28) 



Here R±, Rn, Rq, U are four if-dependent parameter functions. U(K) is a 
4- velocity with only a longitudinal spatial component: 

U(K)= 1 (K)(l,0,0,v(K)), with 1= —L=. (29) 

V 1 — v 



1G 



Its value depends, of course, on t he measurement frame. The "Yano-Koonin 
velocity" v(K) can be calculatedO in an arbitrary reference frame from the 
second central space-time moments of S(x, K). It is, to a good approximation, 
the longitudinal velocity of the fluid element from which most of the particles 
with momentum K are emitted©!]] For sources with boost-invariant longi- 
tudinal expansion velocity the YK-rapidity associated with v(K) is linearly 
related to the pair rapidity 

The other three YKP parameters do not depend on the longitudinal ve- 
locity of the observer. (This distinguishes the YKP form ( |28| ) from the Pratt- 
Bertsch parametrizationHH'ID] which results from eliminating q° in (|27j).) Their 
physical interpretation is easiest in terms of coordinates measured in the frame 
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where v(K) vanishes. There they are given byL 

R\{K) = <y 2 ), (30) 
R\{K) = (Cz-(p l /f3 ± )S;) 2 )-(f3 l /p ± ) 2 (f)^(S 2 ), (31) 

Rl(K) = ({t-5;/l3 ± ) 2 )-(f)/Pi*(P). (32) 

R±, i?|| and Rq thus measure, approximately, the (-fT-dependent) transverse, 
longitudinal and temporal regions of homogeneity of the source in the local co- 
moving frame of the emitter. The approximation in ( |3l| , [32"| ) consists of dropping 
terms which for the model discussed below vanish in the absencej^f transverse 
flow and were found to be small even for finite transverse flowHJ Note that 
it leads to a complete separation of the spatial and temporal aspects of the 
source. This separation is spoiled by sources with (x 2 ) 7^ (y 2 )- For our source 
this happens for non-zero transverse (in particular for large) transverse flow 
r/f, but for opaque sources where particle emission is surface dominated^ this 
occurs even without transverse flowP~ 



4 A model for a finite expanding source 



For our quantitativestudies we used the following model for an expanding 
thermalized source: 113 



S( X ,K)= M ^- Y h, V 



K-u{x) (r-r ) 2 



T(x) 2(At) 



r 

2~R2 



(v-vo) 2 

2(Ar]) 2 



(33) 



8tt 4 At 

1 the spacetime rapidity 77 = \ ln[(f + z)/(t — z)), and the 



Here r — x + y 

longitudinal proper time r = \/t 2 — z 2 parametrize the spacetime coordinates 
x M , with measure d 4 x = rdrdrjrdr d<f). Y = h M(-Eic + Kl)/(Ek — Kl)] and 
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Mj_ = \J m 2 + K 2 ^ parametrize the longitudinal and transverse components of 
the pair momentum K. \/2R is the transverse geometric (Gaussian) radius 
of the source, tq its average freeze-out proper time, At the mean proper time 
duration of particle emission, and A77 parametrizes the finite longitudinal ex- 
tension of the source. T(x) is the freeze-out temperature; if you don't like the 
idea of thermalization in heavy ion collisions, you can think of it as a param- 
eter that describes the random distribution of the particle momenta at each 
space-time point around their average value. The latter is parametrized by a 
collective flow velocity u^(x) in the form 

^(x) = (cosh 77 cosh rit(r), simi^r) e r , sinh 77 cosher)) , (34) 

with a boost-invariant longitudinal flow rapidity rji = rj (y\ = z/t) and a linear 
transverse flow rapidity profile 

m(r) = Vf (£) • (35) 

r)f scales the strength of the transverse flow. The exponent of the Boltzmann 
factor in ( |33| ) can then be written as 

K ■ u(x) = Mj_ cosh(Y — 77) cosher) — K ±-e r s'mhr]t(r) . (36) 

For vanishing transverse flow (77/ = Q) the source depends only on Mj_, and re- 
mains azimuthally symmetric for all Kj_ . Since in the absence of transverse flow 
the /3-dependent terms in ( |3l| ) and (|3|) vanish and the source itself depends 
only on M±, all three YKP radius parameters then show perfect M^-scaling. 
Plotted as functions of Mx, they coincide for pion and kaon pairs (see Fig. El 
left column). For non-zero transverse flow (right column) this Mj_-scaling is 
broken by two effects: (1) The thermal exponent (^) receives an additional 
contribution proportional to K± = v/ M 2 — w? . (2) The terms which were ne- 
glected in the second equalities of ([H],[32|) are non-zero, and they also depend 
on P± = K±/Ek. Both effects induce an explicit rest mass dependence and 
destroy the Mj_-scaling of the YKP size parameters. 

5 iST-dependence of YKP parameters and collective flow 

Collective expansion induces correlations between coordinates and momenta 
in the source, and these result in a dependence of the HBT parameters on the 
pair momentum K. At each point in the source the local velocity distribution 
is centered around the average fluid velocity; two points whose fluid elements 
move rapidly relative to each other are thus unlikely to contribute particles with 
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small relative momenta. Essentially only such regions in the source contribute 
to the correlation function whose fluid elements move with velocities close to 
the velocity of the observed particle pair. 

5.1 The Yano-Koonin velocity and longitudinal flow 

Fig. |l| shows (for pion pairs) the dependence of the YK velocity on the pair 
momentum K . In Fig. |l|a we show the YK rapidity Y YK = | ln[(l +v)/(l — v)] 
as a function of the pair rapidity Y (both relative to the CMS) for different 
values of K±, in Fig. |l|b the same quantity as a function of K± for different 
Y. Solid lines are without transverse flow, dashed lines are for 77/ = 0.6. For 
large K± pairs, the YK rest frame approaches the LCMS (which moves with 
the pair rapidity Y) ; in this limit all pairs are thus emitted from a small region 
in the source which moves with the same longitudinal velocity as the pair. For 
small K± the YK frame is considerably slower than the LCMS; this is due to 
the thermal smearing^ of the particle velocities in our source around the local 
fluid velocity 1^(2:) 113 The linear relationship between the rapidity Y YK of the 
Yano-Koonin frame and the pion pair rapidity Y is a direct reflection of the 
boost-invariant longitudinal expansion flow[!3 For a non-expanding source Y YK 
would be independent of Y. Additional transverse flow is seen to have nearly no 
effect. The dependence of the YK velocity on the pair rapidity thus measures 
directly the longitudinal expansion of the source and cleanly separates it from 
its transverse dynamics. 




Figure 1: (a) The Yano-Koonin rapidity for pion pairs, as a function of the pair c.m. rapidity 
Y, for various values of K± and two values for the transverse flow rjf. (b) The same, but 
plotted against K± for various values of Y and r)f. Source parameters: T = 140 MeV, 
Ar) = 1.2, R = 3 fm, r = 3 fm/c, At = 1 fm/c. 

The NA49 data for 160 A GeV Pb+Pb collisions!!! show very clearly such 
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a more or less linear rise of the Yano-Koonin source rapidity with the rapidity 
of the pion pair. Thi s confirms, in the most transparent way imaginable, their 



earlier conclusion!.^ based on the Y-dependence of the longitudinal radius pa- 



rameter Ri in the Pratt-Bertsch parametrization that the source created in 200 
A GeV S+A collisions expands longitudinally in a nearly boost-invariant way. 

Note that this longitudinal flow need not be of hydrodynamical (pressure 
generated) nature. Similar longitudinal position-momentum correlations arise 
in string fragmentation. This should cause a similar linear rise of the YK- 
rapidity with the pair rapidity in jet fragmentation (with the z-axis oriented 
along the jet axis). It would be interesting to confirm this prediction in e + e~ 
or pp collisions. 

5.2 M±- dependence of YKP radii; transverse flow 

If the source expands rapidly and features large velocity gradients, the "regions 
of homogeneity" contributing to the correlation function will be small. Their 
size will be inversely related to the velocity gradients, scaled by a "thermal 
smearing factor" y/T/M±_ which characterizes the width of the Boltzmann 
distribution^ If one evaluates the expectation values ( [30f]32| ) by saddle point 
integration one finds for pairs with Y = 



R\ = Rl, R 2 = (AU)\ Rj = Ll, (37) 



with 



1 1 1 



(38) 



(At,) 2 = (Ar) 2 + 2(Jri+Ll-r ) , (39) 



2 



1 



LI (t A V ) 2 L\ 



(40) 



ow 



where Rfiow and Lfl ow are the transverse and longitudinal "dynamical lengths 
of homogeneity" due to the expansion velocity gradients: 



^ (M ^ = vj^d^Wr^^ (41) 



L ^ u, - ] : -^ = ^Vi£> (42) 



where ui is the longitudinal 4-velocity. 
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Thus, for expanding sources, the HBT radius parameters are generically 
decreasing f unctions of M±. The slope of this decrease grows with the expan- 
sion ratellSO (this cannot be seen in the saddle point approximated expres- 
sions above). Longitudinal expansion affects mostly the longitudinal radius 
parameter Ru and the temporal parameter i?oF— I the latter is a secondary ef- 
fect since particles from different points arc usually emitted at different times, 
and a decreasing longitudinal homogeneity length thus also leads to a reduced 
effective duration of particle emission (see lower panels in Fig. ^|). 



pion pairs 
kaon pairs 




■■ 0.6 




m^oo 



400 m K 600 



oo m K 600 800 1000 m n 20i 

M ± (MeV/c 2 ) M ± (MeV/c 2 ) 

Figure 2: The YKP radii Rj_, R\\, and Ro (top to bottom) for zero transverse flow (left 
column) and for r/f = 0.6 (right column), as functions of M± for pairs at Y cm = 0. Solid 
(dashed) lines are for pions (kaons). The breaking of the M^-scaling by transverse flow is 
obvious in the right column. For nonzero transverse flow Rq also does not agree exactly with 

the effective source lifetime yj (t 2 ) (lower right panel). Source parameters as in Fig. |l|. 
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The transverse radius parameter R±_ is invariant under longitudinal boosts 
and thus not affected at all by longitudinal expansion (upper left panel in 
Fig. ||). It begins, however, to drop as a function of M± if the source expands 
in the transverse directions (upper right panel). Comparing the lower two 
left and right panels in Fig. g one sees that the sensitivity of Ru and Rq to 
transverse flow is much weaker0 Transverse (longitudinal) flow thus mostly 
affects the transverse (longitudinal) regions of homogeneity. 

While longitudinal "flow" is not necessarily a signature for nuclear collec- 
tivity but could be "faked" as discussed at the end of the previous subsection, 
transverse flow is much more generic in this respect: the ingoing channel has 
no transverse collective motion, and the only mechanism imaginable for the 
creation of transverse collective dynamics during the collision is multiple (re- 
)scattering among the secondaries, leading ultimately to hydrodynamic trans- 
verse flow. 

Unfortunately, the observation of an Mj_-dependence of R± by itself is 
not sufficient to prove the existence of radial transverse flow. It can also be 
created bjj_Qther types of transverse gradients, e.g. a transverse temperature 
gradient oHJ't^l To exclude such a possibility one must check the M^-scaling 
of the YKP radii, i.e. the independence of the functions Ri(M±) (i =_L, ||,0) 
of the particle rest mass (which is not broken by temperature gradients) . Since 
different particle species are affected differently by resonance decays, such a 
check further requires the elimination of resonance effects. 



6 Resonance decays 

Resonance decays contribute additional pions at low M± ; these pions originate 
from a larger region than the direct ones, due to resonance propagation before 
decay. They cause an M j -dep endent modification of the HBT radii. 

Quantitative studies &3c3 have shown that the resonances can be subdi- 
vided into three classes with different characteristic effects on the correlator: 

(i) Short-lived resonances with lifetimes up to a few fm/c do not propagate 
far outside the region of thermal emission and thus affect R± only marginally. 
They contribute to Rq and up to about 1 fm via their lifetime; i?y is larger 
if pion emission occurs later because for approximately boost-invariant expan- 
sion the longitudinal velocity gradient decreases as a function of time. 

(ii) Long-lived resonances with lifetimes of more than several hundred fm/c 
do not contribute to the measured correlation and thus only reduce the cor- 
relation strength (the intercept at q = 0), without changing the shape of the 
correlator^] Decaying at large distances from their production point, they sim- 

a A reduced correlation strength in the two-particle sector could also arise from partial phase 
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ulatc a very large source which contributes to the correlation signal only for 
unmeasurably small relative momenta. 

(hi) Only the u> meson with its lifetime of 23.4 fm/c does not fall in either 
of these two classes and can thus distort the form of the correlation function. 
It contributes a second bump at small q to the correlator, giving it a non- 
Gaussian shape and thus complicatingEJa the extraction of HBT radii from a 
Gaussian fit. At small M± up to 10% of the pions can come from cu decays, 
and this fraction doubles effectively in the correlator since the other pion can 
be a direct one; thus the effect is not always negligible. 
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Figure 3: The inverted q-variance R± and the kurtosis Ax (the index s in the figure stands 
for "sideward") at Y = as functions of K±. Left: r)j = (no transverse flow). Right: 
r\j = 0.3. The difference between dashed and solid lines is entirely dominated by ui decays. 
Source parameters as in Fig. hi except for R = 5 fm. 

In a detailed model studylH we showed that resonance contributions can 
be identified through the non-Gaussian features in the correlator induced by 
the tails in the emission function resulting from resonance decays. To this end 
one computes the second and fourth order g-moments of the correlatorliil The 
second order moments define the HBT radii, while the kurtosis (the normalized 
fourth order moments) provide a lowest order measure for the deviations from 
a Gaussian shape. We found lil that, at least for the model (|33|), a positive 
kurtosis can always be associated with resonance decay contributions (Fig. [| 
left panel). Strong flow also generates a non-zero, but small and apparently 
always negative kurtosis (Fig. ||, right panel). Any Mj_-dependence of R± 
which is associated with a positive M^-dependent kurtosis must therefore be 



coherence in the source By comparing two- and three-particle correlations, the intercept- 
reducing effects of resonances can be eliminated, And the degree of coherence resp. chaoticity 
in the source can be unambiguously determinedQ 
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regarded with suspicion; an Mj_-dependence of R± with a vanishing or negative 
kurtosis, however, cannot be blamed on resonance decays. 

In our model, the first situation is realized for a source without transverse 
expansion (left panel of Fig. 1): At small M± the u) contribution increases 
R± by up to 0.5 fm while for M± > 600 MeV it dies out. The effect on Rj_ 
is small because the heavy u moves slowly and doesn't travel very far before 
decaying. The resonance contribution is clearly visible in the positive kurtosis 
(lower curve). For non-zero transverse flow (right panel) there is no resonance 
contribution to R± ; this is because for finite flow the effective source size for the 
heavier u> is smaller than for the direct pions, and the w-decay pions thus always 
remain buried under the much more abundant direct ones. Correspondingly 
the kurtosis essentially vanishes; in fact, it is slightly negative, due to the weak 
non-Gaussian features induced by the transverse flow. 

7 Opaque sources 

The emission function ( p3| ) is only one of an infinity of possible source para- 
metrizations. Its form permits, however, easy implementation of most of the 
important features of the sources created in heavy ion collisions. Still, there is 
one important physical situation which cannot be parametrized by the formula 
(|33|): if the source emits particles not from the entire volume, but only from 
a thin surface layer. This is how the sun radiates photons, and this is also 
an often suggested picture for the slow hadronization of long-lived QGP blobs 
through a deflagration-type strong first order transition. 

The significance of such a phenomenon for HBT interferometry was real- 
ized by Heiselberg and Vischerl§3 who pointed out that an effective emission 
region which is part of a thin surface layer has a much smaller extension in 
the "outward" or x-direction than in the "sideward" or y-direction. In other 
words, such "opaque sources" have (x 2 — y 2 ) < 0. Depending on the degree 
of opacity (the thickness of the surface layer relative to the source radius) this 
difference can be large and negative. The authors pointed out E0] that this 
leads to the possibility of a smaller "outward" than "sideward" HBT radius 
parameter in the Pratt-Bertsch parametrization, even at K± = 0. Recently 
B. Tomasik showed that in the YKP parametrization opacity effects would 
show up even more spectacularly by a "lifetime parameter" Rq which would 
diverge to — oo in the limit Kj_ — > resp. 13 \ — > (see Eq. (p2])). 

The source ([33]) can be made opaquet3t§ by multiplying it by the factor 
exp (— Kl e s/X)) where A is the mean free path and nl c g is the effective travelling 
distance of the emitted particle through matter in the sourcelH 

Fig. 0c shows the "temporal" YKP radius parameter Rq as a function of 
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M± for sources with different degrees of opacity. For different opacities R/X, 
the transverse source parameters R and rjf were readjustedEl to give the same 
measured single particle slope. The crucial features of opacity are clearly 
visible: the negative contribution ~ (x 2 — y 2 ) in (|3|) drives R 2 to negative 
values at small Kj_, and this happens the sooner the shorter the mean free 
path A, i.e. the thinner the surface layer is. 

Fig. [I|c implies that thin surfaces with A < R are essentially excluded. Pion 
freeze-out occurs not in the form of surface emission, but happens in bulk. 160 
A GeV Pb+Pb collisions are thus quite similar to the Early Universe: an early 
stage of complete opaqueness (here lasting for about 8 fm/c) is followed by a 
rather sudden (At ~l-2 fm/c) transition to complete transparency. In both 
cases this transition is due to the expansion and cooling of the system, causing 
a rapid increase of the mean free path of the particles. In the Early Universe 
photon decoupling is triggered by the recombination of electrons and ions into 
neutral atoms; in heavy ion collisions at SPS energies pion decoupling is caused 
mostly by the rapid cooling and dilution of the baryon density, since baryonic 
resonances with their strong coupling to the pion channel provide the "glue" 
needed for keeping the system close to local thermal equilibrium. 
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Figure 4: The YKP radius parameters Kj_,-R|i, and Rq as functions of Mj_ for pion pairs 

from 158 A GeV/c Pb+Pb collisiontll at slightly forward pair rapidity Y CM = 1.25. Solid, 
dashed and dotted curves correspond to R/X = 0,1,10, respectively. (The corresponding 
pairs (r)f,R) are (0.4,5.74 fm), (0.345,4.83 fm), and (0.215,3.35 fm), respectivelyE3) Other 
source parameters: T=120 MeV, tq=8 fm/c, At=2 fm/c, Ar;=1.3. 




8 Analysis of Pb+Pb data 



Fig. ^ shows a numerical fit 
expressions (|30"|)-(|32|) with our model source (33) 
collaboration in 158 A GeV/c Pb 



of the YKP radius parameters, using the 
to data collected by the NA49 

3£ 



-Pb collisions!^ The fit includes resonance 
decay contributions to the single particle spectra, but not to the 2-particle 
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correlations. Based on sample calculations the latter are, however, expected 
to be inside the systematic error of the data. 

The width of the pion rapidity distribution is reproduced with A?/=1.3. 
To and A?y are then fixedEl by the magnitude of Ru and Rq. The magnitude 
of R±(K± = 0) fixes the radius R once T and r]f are known. The latter are 
obtained from the Mj_-dependence of R j_ , albeit not independently: essentially 
only the combination 77/ a/ Mj_/T (the velocity gradient divided by the thermal 
smearing factor) can be extracted — H This is similar to the single particle 



spectra whose Mj_-slopes determine only an effective blushiftcd temperaturepl 
T c ff = T^J ^ V J . The correlations between T and rjf are, however, exactly 

opposite in the two cases: for a fixed spectral slope T must be decreased if 
rif increases while a fixed M_L-slope of R± requires decreasing values of T]f 
if T is reducedEl The combination of single-particle spectra and two-particle 
correlation thus allows for a separate determination ofT and rjf. 

The fit in Fig. || corresponds to an average transverse flow velocity Vf=Q.4A, 
combined with a freeze-out temperature of about 120 MeV. Similar values were 
advocated by KampferE^j based on a simultaneous analysis of single particle 
spectra of various different hadron species from Pb+Pb collisions. 

Let us discuss in more detail the numbers resulting from this fit. First, the 
transverse size parameter R ss 6 fm is surpri singly large. At mid-rapidity even 



larger values (R ~ 7 fm) are requiredr^fH Resonance contributions are not 
expected to reduce it by more than 0.5 fmEU The transverse flow correction to 
R± is appreciable, resulting in a visible transverse homogeneity length of only 
about 5 fm at small K±, but even this number is large. R—6-7 fm corresponds 
to an r.m.s. radius r rms = yj (x 2 + y 2 ) « 8.5-10 fm of the pion source, to be 
compared with an r.m.s. radius rf^ s — 1.2 x A 1 / 3 * \/2/5 fm = 4.5 fm for 
the density distribution of the original Pb nucleus projected on the transverse 
plane. This implies a transverse expansion of the reaction zone by a linear 
factor ~ 2. That we also find a large transverse flow velocity renders the picture 
consistent. The longitudinal size of the collision region at the point where the 
pressure in the system began to drive the transverse expansion can be estimated 



2\ 



' initial 



as follows: for the source to expand in, say, the y-direction from yjy 
1.2 A x l z j\fh fm = 3.2 fm to \J (j/ 2 ) final = R ~ 6fm with an average transverse 
flow velocity of at most v = 0.44 c (the freeze-out value determined from the fit 
shown in Fig. ||) requires a time of at least (6 — 3.2)/0.44 fm/c ~ 6.5 fm/c. Due 
to the selfsimilarity of the longitudinal expansion the longitudinal dimension of 
the source grows linearly with r. If the total expansion time until freeze-out is 
given by the fit parameter To = 8 fm/c, the source expanded in the ~ 6.5 fm/c 
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during which there was transverse expansion by a factor 8/(8-6.5) = 8/1.5 ~ 5 
in the longitudinal direction. We conclude that the fireball volume must have 
expanded by a factor 5 * 2 2 w 20 between the onset of transverse expansion 
and freeze-out! This is the clearest evidence for strong collective dynamical 
behaviour in ultra-relativistic heavy- ion collisions so far. 

The local comoving energy density at freeze-out can be estimated from the 
fitted values for T and rjf. The thermal energy density of a hadron resonance 
gas at T = 120 MeV and moderate baryon chemical potential is of the order 
of 100 MeV/fm 3 . The large average transverse flow velocity of (vf) sw 0.44 
implies that about 25% flow energy must be added in the lab frame. This 
results in an estimate of about 0.1 GeV/fm 3 x 1.25 x 20 w 2.5 GeV/fm 3 for 
the energy density of the reaction zone at the onset of transverse expansion. 
This is well above the critical energy density < 0.9 GeV/fm 3 predicted by 
lattice QCD for deconfined quark-gluon mattered Whether this energy density 
was fully thermalized is, of course, a different question. It must, however, have 
been accompanied by transverse pressure (i.e. some degree of equilibration 
of momenta must have occurred already before this point), because otherwise 
transverse expansion could not have been initiated. 

9 Conclusions 

I hope to have shown that 

• two-particle correlation functions from heavy-ion collisions provide valu- 
able information both on the geometry and the dynamical state of the 
reaction zone at freeze-out; 

• a comprehensive and simultaneous analysis of single-particle spectra and 
two-particle correlations, with the help of models which provide a real- 
istic parametrization of the emission function, allows for an essentially 
complete reconstruction of the final state of the reaction zone, which can 
serve as a reliable basis for theoretical back-extrapolations towards the 
interesting hot and dense early stages of the collision; 

• simple and conservative estimates, based on the crucial new informa- 
tion from HBT measurements on the large transverse size of the source 
at freeze-out and using only energy conservation, lead to the conclusion 
that in Pb+Pb collisions at CERN, before the onset of transverse ex- 
pansion, the energy density exceeded comfortably the critical value for 
the formation of a color deconfined state of quarks and gluons. There 
is, however, no evidence for long time delays due to hadronization of the 
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QGP, and pions freeze out in bulk rather than from the surface of the 
collision fireball. This is in line with lattice results which predict at most 
a weakly first order confinement transition, and with other evidence cB 
for rapid hadronization. 

I would like to thank the organizers of CRIS98 for their hospitality and support. 
The results reported here were obtained in collaboration with D. Anchishkin, 
S. Chapman, P. Scotto, C. Slotta, B. Tomasik, U. Wiedemann, Y.-F. Wu, and 
Q.H. Zhang. Our work was supported by grants from DAAD, DFG, NSFC, 
BMBF, GSI and the Alexander von Humboldt Foundation. 



References 

1. R. Hanbury Brown and R.Q. Twiss, Phil. Mag. 45, 633 (1954); and 
Nature 177, 27 (1956); 178, 1046 and 1447 (1956). 

2. G. Goldhaber, S. Goldhaber, W. Lee, and A. Pais, Phys. Rev. 120, 300 
(1960). 

3. J. Schmidt-Sorensen, this volume, and references therein. 

4. U. Heinz and Q .H. Zhang, Phys Rev. C 56, 426 (1997); H. Heiselberg 



and A. Vischer, |nucl-th/9707036 



5. K.S. Lee and U. Heinz, Z. Phys. C 43, 425 (1989); K.S. Lee, U. Heinz 
and E. Schnedermann, Z. Phys. C 48, 525 (1990); E. Schnedermann 
and U. Heinz, Phys. Rev. Lett. 69, 2908 (1992); E. Schnedermann, J. 
Sollfrank and U. Heinz, Phys. Rev. C 48, 2462 (1993); E. Schnedermann 
and U. Heinz, Phys. Rev. C 50, 1675 (1994); U. Heinz, in Hot Hadronic 
Matter: Theory and Experiment, (J. Letessier et al., eds.), NATO ASI 
Series B346 (Plenum, New York, 1995), p. 413. 

6. J. Stachel et al, [E814 Collaboration], Nucl. Phys. A 566, 183c (1994); 
P. Braun-Munzinger et al, Phys. Lett. B 344, 43 (1995). 

7. See proceedings of Quark Matter 96, Nucl. Phys. A 610 (1996), and of 
Quark Matter 97, Nucl. Phys. A (1998), in press. 

8. M. Herrmann and G.F. Bertsch, Phys. Rev. C 51, 328 (1995). 

9. S. Chapman, P. Scotto and U. Heinz, Phys. Rev. Lett. 74, 4400 (1995); 
and Heavy Ion Physics 1, 1 (1995). 

10. S. Chapman, J.R. Nix and U. Heinz, Phys. Rev. C 52, 2694 (1995). 

11. S. Pratt, Phys. Rev. Lett. 53, 1219 (1984); Phys. Rev. D 33, 1314 
(1986). 

12. A.N. Makhlin and Yu.M. Sinyukov, Z. Phys. C 39, 69 (1988). 

13. S.V. Akkelin and Y.M. Sinyukov, Phys. Lett. B 356, 525 (1995). 

14. T. Csorgo and B. Lorstad, Phys. Rev. C 54, 1396 (1996). 



28 



15. B. Tomasik and U. Heinz, Eu r. Phys. J. C (23. Feb. 1998) (online 
publication) [ |nucl-th/9707001 1. 

16. U.A. Wiedemann, P. Scotto and U. Heinz, Phys. Rev. C 53, 918 (1996). 

17. U. Heinz et al, Phys. Lett. B 382, 181 (1996); Wu Y.-F. et al, Eur. 
Phys. J. C 1, 599 (1998). 

18. B.R. Schlei et al, Phys. Lett. B 293, 275 (1992); J. Bolz et al, ibid. 
300, 404 (1993); Phys. Rev. D 47, 3860 (1993). 

19. U.A. Wiedemann and U. Heinz, Phys. Rev. C 56, R610 (1997); and 
Phys. Rev. C 56, 3265 (1997). 

20. U. Heinz, in Correlations and Clustering Phenomena in Subatomic 
Physics, (M.N. Harakeh et at, eds.), NATO ASI Series B359 (Plenum, 
New York, 1997), p. 137. 

21. M. Gyulassy, S.K. Kauffmann, and L.W. Wilson, Phys. Rev. C 20, 2267 
(1979). 

22. E. Shuryak, Phys. Lett. B 44, 387 (1973); Sov. J. Nucl. Phys. 18, 667 
(1974). 

23. S. Chapman and U. Heinz, Phys. Lett. B 340, 250 (1994). 

24. D. Anchishkin, U. Heinz, and P. Renk, Phys. Rev. C 57, 1428 (1998). 

25. Q.H. Zhang, U.A. Wiedemann, C. Slotta, and U. Heinz, Phys. Lett. B 
407, 33 (1997). 

26. U.A. Wiedemann et al, Phys. Rev. C 56, R614 (1997). 

27. T. Csorgo and J. Zimanyi, Phys. Rev. Lett. 80, 916 (1998); J. Zimanyi 



and T. Csorgo, |hep-ph/9705432 



28. U.A. Wiedemann, |aucl-th/9801009 



29. D. Pelte an d H. Merlitz, this volume, and references therein | aucl- 
th/9806049 



30. HW. Barz, Phys. Rev. C 53, 2536 (1996); S. Pratt and M.B. Tsang, 
Phys. Rev. C 36, 2390 (1987). 

31. D. Anchishkin and U. Heinz, in preparation. 

32. U. Heinz and U.A. Wiedemann, in preparation. 

33. D.A. Brown and P. Danielewicz, Phys. Lett. B 398, 252 (1997); and 
Phys. Rev. C 57, 2474 (1998). 

34. U. Heinz, in: Proceedings of the 5th International Workshop on Relativis- 
tic Aspects of Nuclear Phy sics, (T. Kodama et al, eds.), World Scientific, 
Singapore, 1998, in press |nucl-th/9710065|1 . 

35. J. Ellis, K. Geiger, U. Heinz and U.A. Wiedemann, in preparation. 

36. I.V. Andreev, M. Pliimer, and R.M. Weiner, Int. J. Mod. Phys. A 8, 
4577 (1993). 

37. H. Heiselberg and A.P. Vi scher, Eur. Phys J. C 2, 593 (1998). 

38. B. Tomasik and U. Heinz, |nucl-th/9805016|, subm. to Eur. Phys. J. C . 



29 



39. H. Appelshauser et al. (NA49 Collaboration), Eur. Phys. J. C 2, 661 
(1998); P.G. Jones et al. (NA49 Collaboration), Nucl. Phys. A 610, 
188c (1996); R. Ganz, this volume. 

40. T. Alber et al. (NA35 Collaboration), Z. Phys. C 66, 77 (1995). 

41. U.A. Wiedemann, B. To masik, and U. Hei nz, in Quark Matter 97, Nucl. 
Phys. A (1998), in press |nucl-th/9801017 l. 

42. S. Schonfeld er, PhD thesis, M PI ffir Physik, Miinchen (1996). 

43. B. Kampfer, [nucl-th/9612336| , and private communication. 

44. E. Laermann, Nucl. Phys. A 610, lc (1996). 

45. J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank, and J. Rafelski, Phys. 
Rev. C 51, 3408 (1995). 



30 



